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We discuss in more details the theory of low energy excitations in quantum impurity 
problems with an external field at vanishing temperature, giving further support to results 
of the previous paper. We then extend these results to the next order at low frequency, 
obtaining in particular the exact expression, as a function of the bias, of the first two 
derivatives of the response function x"i^) at a; = in the double well problem of dissipative 
quantum mechanics. 

We also extend our approach to the case of non vanishing temperature and no external 
field. Fendley et al. had obtained in that case an expression for the Hall conductance with 
a single impurity using a Landauer-Biittiker type approach. We recover their result in the 
framework of linear response theory using renormalized form-factors. We also obtain, as 
a function of the temperature, the first derivative of the response function x"{^) at w = 
in the double well problem of dissipative quantum mechanics. 
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1. Introduction 



A lot of progress has been accomplished recently in the non perturbative computation 
of correlation functions for integrable quantum field theories in 1+1 dimensions 
Interesting physical applications have been found in particular in one dimensional quan- 
tum impurity problems [Q, 0,0. Unfortunately, many of these computations have been 
restricted to the case of vanishing temperature and external field, preventing in most cases 
comparison with experimental data; a notable exception concerns the DC properties for 
tunneling between edge states in fractional quantum Hall devices f^. 

The introduction of a temperature or an external field, while it does not break integra- 
bility, renders the computation of correlators using form-factors more difficult (putting a 
temperature seems somewhat more natural in the approach of , which however produces 
only rather implicit results). Few works have addressed this question so far, except in the 
free fermion case, where however the problem is already quite non trivial for spin operators 
1^ . We note that some interesting results have also been obtained in the case of interacting 
massive bulk theories in the limit of very low temperatures. In (henceforth refered to 
as I) we have considered quantum impurity problems at T = but with an external field. 
We were able to determine two important low energy properties, the limit lim^j^o ^ ^ of 
the dynamical susceptibiliy in the double well problem of dissipative quantum mechanics, 
and the component of the noises in a four terminal geometry for tunneling between 
edges in fractional quantum Hall devices. Some global properties, like the existence of 
potential singularities, were also examined. 

The first purpose of this sequel is to place the low energy analysis of I on firmer 
grounds, and to compute the next order at low frequency, eg the uJ^ term in ■}(^' {lS). This 
is the subject of sections 2 and 3. 

Our second purpose is to carry out a similar analysis for the case of a non vanishing 
temperature (section 4). We discuss the formula for the conductance through a point 



contact that had been obtained previously by Fendley et al. using a Landauer-Biittiker [|10 
type approach, and recover it, at the price of some reasonable hypothesis, in the framework 
of linear response theory. Building on this analysis, we obtain the limit lim^^^o ^ of 
the dynamical susceptibiliy in the double well problem of dissipative quantum mechanics 
at finite temperature. 
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2. A closer look at low energy excitations with T = 0,V ^ 

This paper is a sequel to I and the notation introduced there wiU be used here. The 
hamiltonian we study, 



consists of a free bulk theory with a boundary interaction. The detailed form of this 
interaction, Hb, is different for the dissipative quantum mechanics and the tunneling 
problems. 

The principle of our approach using massless scattering is described in full details 
in I. In brief, the bulk degres of freedom are described by a massless sine-Gordon theory 
which has solitons/anti-solitons and bound states as fundamental "bare" excitations. In 
the presence of a voltage, the ground state consists of a Fermi sea filled with solitons, in a 
way easily controlled using the Bethe ansatz. Physical properties of interest are determined 
by the structure of excitations over this voltage dependent ground state. The boundary 
interaction is finally taken into account by the introduction of refiection matrices. 

In this section we discuss the voltage dependent ground state and its excitations more 
thoroughly. The main problem we are addressing is the scattering of these excitations, 
which is "dressed" by the presence of the Fermi sea, and the resulting dressed form-factors. 
Most of the discussion is centered on the scattering at coincident rapidities, where, for non- 
neutral excitations, a non trivial phase seems to appear. This phase, and its possible effects 
on the low energy properties, are discussed in details. 

2.1. The Fermi sea 

Recall that at zero temperature and in the presence of an external field, the ground 
state contains only solitons. The model is massless and a convenient way to use the Bethe 
ansatz in this problem is to "unfold" the half-line problem to obtain a problem on the full 
line. Then we can, for instance, consider a system with only right movers. The system is 
therefore described by the Bethe ansatz equations (setting ^ = 1) 




(2.1) 




(2.2) 
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where we have reinstated an arbitrary energy scale /U (set equal to one in I) and 5 is the 
shift coming from the soliton-soliton S matrix. In the ground state, solitons fill the interval 
6 E [—00, A]0, with a density p obeying 



pio) = ^+ r He~e')p{e')de', (2.3) 

where = S = e^*'^'^, $(6') = ^ • It should be pointed out that the quantization 
condition ( p.2| ) should also involve the reflection matrix but its contribution to the density 
is of order 1/L, and this is negligible for our purposes. Equations of the type ( pT^ ) for a 
general function / 

f{e) = g{e) + [ - e')f{e')de' 

will occur repeatedly in the subsequent analysis. Let us introduce the functional K such 
that ^ 

k*f{e)= I ^{e-e')f{e')de' 

J — oo 

and similarly 

i.f{e) = f{e) 

Then the previous equation reads {I — K) • f = g, and the solution follows by introducing 
an operator L such that 

{i + L){i-k) = i, 



that is 



f{o) = g{e)+ r L{e,9')g{e')de'. 

J — oo 



The functions L is symmetric of its two arguments L(^, 9') = L{9' ^ 9) and can be written 
formally as 

L{9,9')=^{9-9')+ f ^9 - 9")^{9" - 9')d9" + . . . . (2.4) 

J — oo 

It is important to stress that L is not a function of the difference of its arguments, as the 
previous equation clearly shows. Observe that one has (/ + L)K = L, or more explicitly 

/ [6(9 - 9') + L{9, 9')]^{9' - 9")d9' = L{9, 9"). 

J — oo 



^ As in I we use rapidities to parametrize the energy and momentum of a particle. 
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In the particular case of the density we have p{e) = ^ + L{e, 9') e^' dO' . For 
9 < A, p is the density of sohtons in the sea. For 9 > A, p is the density of holes of 
solitons above the sea. The function p can be written maybe more explicitely as an infinite 
series of exponentials using Wiener Hopf integration techniques, see I. We recall the value 
p{A) = ^^J2g. The Fermi rapidity A is determined self consistently by introducing the 
quantity 

e{9) =pe^ $(6' - 9')e{9')d9', (2.5) 

and requiring e{A) = 0. Observe the identity 

de{9) 



d9 

2.2. Neutral excitations 



2ttp{9). (2.6) 



Let us now discuss excitations over the ground state. First, consider a particle hole 
excitation. In order not to make confusion with the other types of particles (eg antisolitons) 
which also appear in this problem at high energy, we call a soliton above the sea a volton, 
and a hole in the sea an antivolton. Suppose the volton has rapidity 9p and the antivolton 
rapidity 9h. The particles are interacting and this induces a shift of the rapidities in the 
seal a rapidity equal to 9a initially becomes 9a + S^'^^9a with the conditions 

ZTT ^ — ' 



P 

=t^e'-+'''''- + y 5{9a -9p + d^'^9a - d^'^9s) (^.7) 
+ Si9a-9p)-S{9a-9h). 



We then define the shift function, describing the change in the Fermi sea, by Lp{9)5^^^9 = 
F{9\9p, 9h). By standard manipulations [0, one finds the equation obeyed by the shift 

F{9\9p,9h)- r ^9-9')F{9'\9p,9^,)d9' = 6{9-9^,)-6{9-9p). (2.8) 

J —oo 

A formal solution of this equation follows as 

F{9\9p,9h)= ['"1(9,9^9'. (2.9) 



The change of energy of the system with this excitation is given by the bare energies of 
the particle and the hole and the contribution from the shift of the Fermi sea 

5(2)^ = ^jie^^ - fie^'^ + / F{e\ep, Oh) ^e^dO, (2.10) 

J — oo 

and this can be shown to coincide with 

S^^^E = e{ep)-e{eh), (2.11) 

with e defined in ( |2.5| ). One has e{9) < for 6* < A, so the energy of hole excitations 
is positive as it should. We note the value e(— oo) = —ngV. Similar computations can 
be carried out for the momentum, with the result that, for these neutral excitations, 
= hence the excitations are massless, as physically expected. 

2.3. Non neutral excitations 

Although physical excitations will always be made of volton-antivolton pairs, it is 
necessary, to understand the scattering properties, to consider excitations with only one 
volton or one antivolton. Such excitations will involve a "Fermi" momentum, that is 
= 5^^^ E±pf. Assuming for the moment that the excitation energy of a volton (resp. 
an antivolton) is still given by e (resp. -e), let us now determine the Fermi momentum. 
Recall that the proper definition of p is via the phase shift collected when a particle goes 
around the world. For a volton excitation this phase reads 

Lp = Lfie^ + 27r ^ 6{9 - 9^^^) + 27r ^ 6{oo). 

a a 

Here, we have taken the usual phase shift 5 = t^— In 5" where 5" is the soliton soliton sine- 
Gordon S-matrix. However, when taking the massless limit of the sine-Gordon model, 
left and right movers do not become totally independent. There remains a RL constant 
scattering phase, 6{oo). This phase appears very rarely in computations, and its meaning is 
not totally clear. But when we pass a particle through the system, this RL shift contributes 
to the momentum, and this is the meaning of the last term in the foregoing equation. At 
leading order, we can neglect the shift of rapidities in the sea and reexpress this as 

Lp = Lue^ + 27r / [5(9 - 9') + 5{oo)]p{9')d9' 

J —oo 



Then, using the relation p = we can write 

Lp = Le{e) + L^ + 2fiLgV5{oo). 

One checks that 5{oo) = —5{—oo) = | ~ ^5 and it foUows that 

Pf = pgV. (2.12) 

Consider now the operator destroying a right moving volton and creating a left moving 
one. Its correlation function, from this extra Fermi momentum, will be alternating with a 
cos2pFX part. On the other hand, the alternating part can be computed by observing that 
the potential V can be absorbed by a redefinition of the field (I)l,r ^ 'Pl,r ± gVx in the 
original hamiltonian. This leads to the identification of this operator with expi{(pL—(pji), of 
conformal weights h = h = g. More generally, the operator expia{(pL — ^Pr) has conformal 
weight h = ga^ and its correlation function involves a part cos2agVx. This corresponds 
to a charge Q = 2ga, in conventions where the soliton has unit charge . 

Notice that the LR phase shift should appear in ( |2.2| ). However, as long as the number 
of L particles is a constant, its presence would simply shift all the by a constant, without 
changing the results for densities or excitation energies. Things are different when the 
number of particles is changed. For instance, adding a R soliton adds up a phase 5(— 00) 
to the rhs of Bethe equations for L movers. Without any correction term, this phase in turn 
will move the Fermi sea, resulting in an induced L charge, an effect which is not expected 
for a purely R excitation. This means that, to observe physical non neutral excitations, 
one has to complement the addition or removal of a particle by an additional phase shift, 
corresponding presumably to changing boundary conditions. Determining the value of this 
phase is not obvious and, in the rest of this section, we embark in a long discussion of the 
properties of the system with this phase to see its effect. Consider for instance a volton 
excitation, associated with an additional phase shift exp2z7r5;,. The following holds 

27r ^ 

/o + 5, =^e'^+^''''^ + T 5{e^ -e, + 5(1)^. - 5('%) (2-13) 
2n ^ 
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with 9p the volton rapidity. By the same logic as before, defining Fb{9\ep) = Lp{e)5^^^e, 
one has ^ 

F^{e\ep) -I $(0 - e')F^{e'\ep) = 5, - 5{e - Op). (2.14) 

J —00 

If we were to add an antivolton, the additional phase shift would be exp — 2indb, so of 
course for pairs voltons antivoltons the total shift would be as before, as one has 

F{e\ep,eh) = FtieiOp) - FtieiOh) 



for any 6b, so ( |2.10| ),( |2.11[ ) hold indenpendently of this phase. Now, for non "neutral" 
excitations we have to require more, namely that the excitation energy also is given by 
5(1)^; = e(^p). Let us see whether this can be satisfied: On the one hand, the solution of 
( pT5| ) reads 



' - 2 + 



ne^L{ep,e)de 



V 



L{9p,9)d9. 



(2.15) 



On the other hand, the excitation energy due to the volton can be written 



fie'- - ^ + / fie''Ft{9'\9p)d9' - ^Fb(-oo|^p). (2.16) 

^ J -oo ^ 



The last term was not present for neutral excitations. First, using the formal solution 
( P^) , one checks that Fh{—oo,9p) is in fact independent of 9p. For neutral excitations, 
this term would be cancelled by its counterpart due to the antivolton. The meaning of 
Fb{—oo,9p) is simple: it is (minus) the number of solitons that are pushed "out" of the 
Fermi sea in the presence of an added soliton above the sea. While these particles have a 
vanishing kinetic energy, they have a non vanishing potential energy, which must be put 
by hand. 

With the shift 5^, the solution of ( |2.14]) reads 



Fb{9\9p) = [5b - 5(00)] 



1 + 



L{9, 9')d9' 



+ 



L{9, 9')d9'. 



[2.11) 



Prom the defining equation for F^ (|2.14D , we can also deduce 

d 



Fb{9\9p) 



d9 

^^Fb{9\9p) =1(9, 9p). 



L{9,9p)-L{9, A)Fb{A\9p 



(2.18) 



Integration by parts in ( |2.16|) givves rise to the alternate expression 

V 



e{ep) = fie'--- + ^ie''Ft{A\e, 



+ 



l-A y 

/ ne^L{e,ep)de + Fi,{A\ep) / ne^L{e,A)de- -Fb{-oo\ep). 

J — CO J — oo 



(2.19) 



Using e{A) = 0, the equaUty e = e wiU hold iff 

-A 



FbiAlOp) 



1 + 



L{A,6)d6 



-Ft{-oo\9p 



A 



L{ep,e)de. (2.20) 



The left hand side is linear in 5h with the coefficient 



1 + 



L{A,e)de 



- 1 



A 



L{-oo,e) de 



(for the second term, the upper integration bound has no importance since the region 
where the integrand is finite is concentrated around — oo). By using the equations for p 
and e together with the values of p[A) and e(— oo) one finds the results 

I L{A,e)de -1 

J —00 



A 



(2.21) 



L{-oo,e)de =2g - 1. 



Hence the left hand side of ( |2.20|) actually is independent of 5^. Let us chose the particular 
value b* of Sb such that F^* (A, A) = 0. Using (|2.18| ) it follows that in that case F^* {9, A) = 
-Fb* {A, 6). One has then 



Fb*{ep\A) 



J^jL{6p,6)-LiA,6)]d6 



1 + j^^ L(A, e)de 



(2.22) 



from which ( g3o|) follows together with Fb* (-oo, 6) = J^^ L{A, 9)d9 = ^ - 1. Hence, 
indeed, the function e gives the excitation energy of non neutral excitations too, and this 
for any value of db- 

Another quantity of interest is the charge of excitations. For a volton for instance, it 
is equal to +1, the bare charge, minus the number of solitons ejected from the sea at — oo, 
which reads —Fb{—oo\9p). So we have 



q±{e) = ±[l + Fb{-oo\e)] 



(2.23) 
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where from now on, we designate voltons and antivoltons by a label ±. The modulus of 
this charge is a constant, so neutral excitations still consist of pairs voltons-antivoltons, as 
physically expected. For a generic choice of shift 5^, one finds that 0016*) = 2g[5b — 
5(cx))] + 2g-l. For this value, one has Fb{A\A) = y/2g[5h - 5(oo)] + ^/2g - 1. 

At this stage, it is useful to consider the effect of shifting the edge of the Fermi sea. 
We would like to compute the corresponding change in energy to order 1/L and see if 
that will furnish more constraints on 5h- So far, we did all computations without worrying 
about such terms. It is for instance not possible to use straightforwardly the formula for 
excitation energies e when one wants such corrections. The safest is to go back to original 
definitions. The energy reads as a discrete sum. When replacing the sum by an integral 
in the ground state, the Euler-Mac Laurin formula can be applied. It shows that there is 
a term proportional to L, no term of order one (as physically expected when there is no 
boundary nor impurity), and a term of order 1/L that determines the central charge. If 
we shift the edge by a small amount 5 A, this term of order 1/L will have a variation of 
order 1/L^ which we do not look for. Only the variation of the extensive term is therefore 
of interest. We have 

where 

oA+SA 



As in P] we can easily rewrite 



E ^ '■^+'^ 



L 2tx 

and a similar equation for non tilde quantities, where 

A+SA 



e{e) =^e^ --+ $(6* - 9')i{9')de' 

^ J -oo 



To get a corection of order 1/L to the energy, we expand E and e as functions of 5 A. One 
finds 

dm Q 



d5A 



d'^6A 



SA=0 

= e(A)L(0,A), 

5A=0 
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from which it follows that the first derivative of E at SA 
derivative one has 



1 d^E 



L d^6A 
It follows finally that 



5A=0 



ZTT 



V v_ 



vanishes, and for the second 



L{9, A)d9 



E-E 



(2.24) 



Now suppose we want to create an excitation of charge Q. For each volton added, a 
certain number of voltons leave the Fermi sea, so for dn voltons added, we have a charge 
Q = [1 + Fi,{—oo\A)]dn. By which amount does the Fermi edge shift if we add one 
volton with minimum possible energy? First, the rapidity of the soliton immediately 
below the added volton shifts, by an amount Si = Second, the added volton must 

lie immediately above, at a rapidity differing from this one by 82 = Xp|^- This leads to 



~ Lp{A) ^" Lp{A) thus a change of energy, using (|2.24| ) 



^ ^ _ 1 gV^ f F,{A\A) + 1 
L 8n \ p{A) 

If we add dn voltons, the change of energy will thus read 

E-E=y[l + Ff,{A\A)]\dnf. 
Using the formulas for the charge, together with the usual formula between the gaps and 



the conformal weights [|TT| leads to 



h 



^[1 + F,{A\A)f = ^[1 + F,{-oo\A)f{dnf = ^. 



(2.25) 



This agrees with what is expected from the Lagrangian. Notice that, once again, the 
result holds for any value of db'- so far, we have obtained no constraint for this unknown 
parameter. Its value is however crucial for the stattering theory, as we now discuss. 



2.4- Renormalized S matrix 

To proceed, we consider the S matrix of excitations. Take for instance two voltons 
at rapidities 6*1, 6*2, with 9i > 62- With the first volton only, the rapidities in the sea 
are shifted to 6^a \ with both particles they are shifted to 6^a'^\ Passing the first volton 
through the system in the presence of the (shifted) Fermi sea only, one gets a phase ^1, 
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while passing it through the system in the presence of the (shifted) Fermi sea and the 
second particle, one gets a phase (^12. These two phases read respectively, 

(2.26) 

4)12 =Le'' + 27r 5^ 5(^? - ^i^^)) + §{61 - 62). 

ct 

Straightforward computations give 

(pi2-(pi = -27:F,{ei\92), (2.27) 

wiht F given in ( |2.17| ). Here, we have assumed that with n voltons, the boundary conditions 
will be changed so that a term ndb will be added to the left hand side of quantization 
equations (similar to (|2.13| )) for solitons in the sea. We then define the S- matrix of voltons 
by 

^++(^i,^2) = -exp[-2z7rF,(^i|^2)], ^1 > ^2- (2.28) 

The previous computation does not require 9i > 62 to be algebraicaly correct. What 
does require this inequality however is the identification of the 5" matrix with the physical 
process of passing the particle 1 through the system. 



We now build a Faddeev Zamolodchikov [0 algebra for our excitations, setting 

1^1,^2 >++= ^++(^1, ^2)1^2,^1 >++, ^1 > 92. (2.29) 
From this, we see that 

1^1,^2 >++= ^;i(^2, ^1)1^2,^1 >++, ^1 < 92. (2.30) 
We will henceforth set 

^++(^1,^2) = -exp[-2z7rFfe(^i|^2)], ^1 > ^2 

(2.31) 

= -exp[2z7rF,(02|^i)], ^2 > ^i- 

The S-matrix satisfies 

^++(^i,^2)5++(^2,^i) = 1. (2.32) 

It is singular at 61 = 6*2, and is not a function of 61 — 62 only. The same analysis for 
antivoltons results in 

5__(^i, ^2) = - exp [-2z7rF6(^i|^2)] , ^1 < ^2 

(2 33) 

-exp[2z7rFb(^i|^2)], ^1 > ^2- 
11 



Observe that the phases are opposite to the ones in ( |2.31| ). This is because the momentum 
of an antivolton, p = Pf + e(^), is a decreasing function of 9 for 6 < A. Finally, one can 
also scatter a volton and an antivolton. One has 



1^1,^2 > + -= (^1,^2)1^2,^1 >- + , (2.34) 

with 

^+_(^i,^2) = -exp[2i7rF6(^i|^2)], (2.35) 

together with 

^_+(0i,02) = -exp[-2z7rF,(^2|^i)]. (2.36) 

In the foregoing equations, there is an additional minus sign compared with the phase shifts 
like ( |2.27| ) because our definition of the S-matrix is through the Faddeev Zamolodchikov 
algebra []T^. Indeed, consider for instance a pair of particles in some scattering theory. 
If (j)i2 is the phase one gets passing the first through the second (this 0i2 should not 
be confused with the one in equation (|2.27|) ), the coordinate wave function reads, before 



symmetrization or antisymmetrization 

gi<^i2/2gZ(pixi+p2X2)^ Xl > X2 

from which 

\p^P2>= j e-^^i2/2e*(Pi^i+P2^2)^+(^^)^+(^2)|0 > 

J Xl <X2 
J Xl >X2 

= -e-"^^'\p2Pi >, 
where we used that the x operators are fermionic, and (j)i2 + 4>2i = 0. 

2.5. Form-factors 

Still assuming a generic Sh-, let us consider the physics at low energy. Right moving, 
low energy excitations are made of voltons antivoltons pairs near the Fermi rapidity 6 = A. 
Because there is an energy scale V , the theory is not relativistically invariant, even though 
neutral excitations have dispersion relation e = p. This means, as observed above, that 
the S-matrix does not depend only on the ratio of energies. However, very close to A, the 
ratio e/V goes to zero while the S-matrix elements go to constants. In that limit, we deal 
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again with a relativistically invariant theory, now depending on Sf,- The key parameter 
now is 

K = Fb{A\A) = y^[5i, - 5ioo)] + ^/2^ - 1. (2.37) 

Since the theory is relativistic, it is convenient to parametrise the energy in terms of a 
renormahsed rapidity (3 (see I). For voltons, we set e(^) = /ie^, 6 > A, and for antivoltons, 
we set e{9) = — /ue^, 9 < A. We normahze asymptotic states to 2tt5{(3i— (32). Note that due 
to the minus sign, for antivoltons, the function P{9) is decreasing so increasing momentum 
means increasing renormahzed rapidities. 

All form factors are now expected to depend on differences of /9's, which are renor- 
mahzed rapidities. The choise 6 = 6* or k = is the only one for which the S-matrix is 
analytic at coincident rapidities. This by itself is probably enough to dictate that choice. 
However, we wish to accumulate more evidence, and keep studying the generic case first. 

We discuss the form factor < 0|j(2)|/32, /3i > \-, where the right hand side physically 

describes a hole particle pair and z = x — t. The z dependence of this correlator follows 
trivially from kinematic considerations, exp niz[e^'^ + e^^], so we restrict to the case ^ = 0. 
We call this form factor /_| (/3i, /32)i. The first axiom to be satisfied is 

/+_(A,/?2)^+-(/3i,/?2) = /-+(/32,/3i) (2.38) 



with the 5" matrix given by its low energy limit described by ( |2.37| ). An obvious solution 
to (|2.38| ) is simply to set 

/+_(/3i,/32) oc exp[-i7rK] 
/_+(/32,/3i) oc -exp[inK\. 



(2.39) 



(2.40) 



The second axiom is 

/+_(/3i, /32 + 2i7r) = ^+_(/3i, /32)/+-(/3i, /32) 
/+_(/3i +2z7r,/32) =S-'_{(3i,(32)f+-{(3u(32). 

And a solution to ( p.40|) is of the form 

/+_[A,/32]oc^e^^/V^/2(e^-/^0^ (2-41) 

where we have put the correct dimension 1/ length of the form factor. Apart from the 
dimensional factor, the form factor is expected to depend on the massless relativistic 



The notation is similar to [§], < 0|i|/32,/3i >= /(/3i,/32). 
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invariant s = e^^ We see from ( 2.41 ) that there is a cut in the s plane along the 



negative real axis. With k 7^ 0, such cut is unavoidable since the S matrix is singular at 
coincident rapidities. 

The last axiom in the relativistic case stems from crossing: 

_ < /3i|i(0)|/32 >-= - z7r,/32). (2.42) 

Leading to the minimal conjecture for the two particle form-factor 

/+_(/3i,/32) = iiiq-e-'^^ef^/^ef"'/^ (e^^"^^)" . (2.43) 

The normalisation of the form factor is related to the charge term which originates from 
crossing, as we now demonstrate. Indeed, we must satisfy the relation 

/oo 
- < PilJizm >- dz = 2n5{(3i-(32)q-. (2.44) 
-00 

After crossing, performing the integral over z, constrains /3i = (32 and brings a jacobian. 
The ratio of e's becomes equal to unity and disappears. The only remaining part, therefore 
is the charge as desired. In the case k = 0, 5^ = 5b* = —\ + — (j_ = —y/2g, so 
one has then 

/+_(/3i,/32) = -z^v^e/3i/2e/3^/2. (2.45) 

Note that instead of (|2.39| ), (|2.38| ) could be solved also by setting, 

U.{(3,, (32) ^ [e^^ - e^^Y ^ ^ 

2.46 

/_+(/32,/3i)oc-(e^^-e^0 • 

However, in the computation of the charge, this would result in an additional factor e'^^ 
in the right hand side of ( p.44|) , which is impossible. 



The same logic can be applied to form-factors with higher number of particles. A 
reasonable guess for the four particle form-factor is 

< 0|j(^)|/34,/33,/32,/3i >— ++0C q-ix{e^'-^ - 1) (e^3+/34-/?i (g/^i + -f e^^ + e^*) 

sinh sinh 



cosh cosh cosh cosh ' 

(2.47) 

Form factors where the /9's have been exchanged will take similar forms, up to phases 



arising from the S-matrix as in (|2.39| ) - reproducing these phases by terms depending on 

14 



the differences of the energies would make the pole axioms and the charge computation 
impossible as in ( |2.46| ) . Note that these phases are singular at coincident rapidities (3i = (32 
or (3s = (34, since the S matrix is. 

Expression ( pj.47D vanishes when two voltons or antivoltons coincide, as physically 
required in the Bethe ansatz. It also has a simple pole as a volton and antivolton annihilate, 
eg (3^ = (32 + iTT, with a residue that is proportional to the two particle form-factor in 
agreement with the axioms in (it is not completely clear how exactly to generalize the 
pole axiom, and therefore what this residue exactly should be). 

Let us now use these form factors to discuss physical quantities, relying heavily on I 
for the definitions of the impurity problems and the quantities under study. The important 
point about ( |2.47|) is that it transforms as FF e'^FF when P ^ P + a. As a result 
it will, together with all the higher form factors, contribute to the noise computation at 
order uj. While the normalisation of successive form- factors are dictated by LSZ reduction 
formulas ]![] , and the normalisation of the first form- factor by the charge, a very non trivial 
sum rule has to be satisfied to ensure that the current correlator reproduces 

S{oo) = J e^-*({jW,j(0)}) = ^, (2.48) 

(see I). We suspect this sum rule would not be satisfied away from n — 0, but we have no 
proof of this. 

Some restrictions are obvious however. The two particle form factor by itself gives to 
the noise a contribution (assuming a; > 0) 

S^'^ = ^{q-r fdx = -^il/2 + g + 2g6,r^. (2.49) 

2n Jq \uj — X J sm2K7r 27r 

The ratio involving k is clearly greater than one. Hence one needs (1/2 + (7 + 2g6h)'^ < 2g. 
For instance, for g < 1/2, this excludes the choice db = 0. The choice 5^ = 1/2 is always 
excluded. In fact, any choice db = est is excluded for g small enough, so db must have some 
non trivial dependence on g. 

2.6. x" revisited. 

Following I, let us now add a boundary interaction. Without a voltage, the effect of 
the boundary is taken into account by refiection matrices R (solutions of the boundary 
Yang-Baxter equation) for the "bare" excitations. With the voltage, one obtains new, 
dressed refiection matrices TZ, for the excitations over the Fermi sea. 
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Consider again the computation of the low frequency behaviour of in dissipative 

quantum mechanics along the same lines as in 1. For a generic value of the two particle 
term would be proportional to, instead of eq. (5.8) in I, 

\R{x)n\x-u)-\\y-^—\ dx. (2.50) 
Change variables x ux to get 



UJ 



1 / \ 2k. 

' X ^ 







[n{ujx)'R*{xu -u)-l] \ J— ^ ) dx. (2.51) 



We are interested in the first non trivial real term, which occurs from the expansion of the 
integrand to order u'^ . Setting IZ = e*'^, the bracket gives a total of four terms 

-^(<P')W[x' + - 1)' - 2x(x - 1)] = 
so the first contribution is equal to 

Observe that this integral is exactly the one that would appear in the computation of the 
current correlator without impurity ( p.49|) : the impurity dependence {4>') and the energy 
dependence actually factor out. This result is true at every order. For instance the next 
order term would we proportional to 



puj puj puj 

/ / / \f{xi,X2,X3,UJ-Xi-X2-Xs\ 

Jo Jo Jo 



[n{xi)n{x2)n*{-xs)n*{xi + X2 + xs-u;)-l] dxidx2dxs 



XiX2Xs{(^ — Xi — X2 — Xs) ' 

where / is the four particle form- factor (|2.47 ). Now the same integral without the term in 
brackets is the one appearing in the computation of the noise, the fourth order term of a 
whole series that sums up to giv. Contributing to the linear term of x" we have to extract 
the term of third order in uj. Make a change of variables x ujx. Using homogeneity of 
the current form- factors, we can rewrite the integral as 

^ |/(xi,a;2,X3, 1 - xi - X2 - xal^ 

Jo Jo Jo 

{lZ[uJXi)lZ[uJX2)l<- {—UJXsj/Z [uJ[Xi -f X2 + X3 — IJ — 1} 



XiX2X3(l — Xi — X2 — Xs) 
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To get x", we expand the integrand to order w^. Setting again IZ = e*"^, the corresponding 
terms are identified, setting X4 = 1 — xi — X2 — X3, as 



Hence the term (</>') ^u;^ factors out, and the remaining integral is again the same as 
the one arising in the computation of the correlator without impurity. Since the form- 
factors are normalized to sum up these integrals to the same term quj/ti for any 5b, we get 
the same formula for the limit (as a; 0) of ■x"{uj)/uj for any value of 5^! Observing that 



^ 4s ln7?,(6')| „ . , where a = -Tja\ is given in I, we recover 



X"{uj) ( -0 d 



lim = e-' — ln7^(6') 

c^^o uj 4gn^ \ dO 



(2.52) 



This is true for any value of db- 
2.7. Summary 

The analysis of the theory of excitations over the new ground state, when done care- 
fully, turns out to involve an extra phase Sb, which has a key infiuence on the dressed 
scattering theory, and thus on the response functions in the presence of a voltage. In this 
long analysis we have looked at various physical quantities to check whether we could fix 
6b. The computation of lim^^^o ^ ^'^ and the corresponding Shiba relation of I are not 
enough to settle the value of db'. in fact, any value would give the right result, although 
the value that leads to k = is the most reasonable physically, and actually the only one 
for which the sum rules for the charge and noise can be explicitely checked. In contrast 
with dissipative quantum mechanics, the value of k is crucial for tunneling between edges 
in quantum Hall devices, also discussed in I. Indeed, if k was non zero and the theory 
admitted multiple particle, low energy excitations, the noise in the presence of an impurity 
would involve terms of an arbitrarily high degree in the refiection matrix. The same fea- 
ture would be observed in the conductance, as will be discussed later in this paper. This 
is in contradiction with formulas that have been derived using a variety of approaches (in 
particular a Boltzmann type equation [||], or a Landauer-Biittiker approach [0, [0), and 



checked against duality, perturbative expansions, numerical and experimental data. To 
summarize: while we have no proof that k = 0, there is a large body of practical and 
conceptual evidence that it is so. 
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3. Beyond the low energy behaviour 

Let us thus assume from now on that 5h = Sh* has been chosen so that k = 0, and that, 
at low energy, the theory is made of free fermions. What can we say for higher energies? 
Let us restrict to energies w < gV so the only processes involved concern volton-antivolton 
pairs. Since excitations have finite energy, we turn back to the original variables 9. We 

consider the form-factor < 0\j{z)\92, Oi > |_. The z dependence is now exp iz[e{9i)—e{92)]- 

Since the S- matrix describes exchange properties, it is natural to expect that the analog 
of (|2.38|) will be satisfied (we denote Fi,* = F), 



/+_(^l,^2)^+-(^l,^2) =/-+(^2,^l). (3.1) 

Generalizing ( |2.39| ), this has the simple solution 

/+_(^i,^2) ocexpH7rF(^i,^2)] 

(3.2) 

/_+(^2,^i)oc-exp[z7rF(^i,^2)]. 

The analog of crossing is also clear. While voltons are originally defined for 9 > A and 
antivoltons for 9 < A, one can continue the associated amplitudes by analyticity through 
the threshold A. Creating an antivolton at 9 > A will simply mean, as before, destroying 
a soliton (called a volton here) at this rapidity. In other words, we expect the analog of 
(S): 

- < ^i|j(0)|^2 >-=< 0b-(0)|^2,^i >-+= f+-{9u92). (3.3) 

The shift of in was necessary in the (3 parametrization of the relativistic case to describe 
both sides of the thresholds, something we simply accomplish here by continuation through 
the threshold. 

Away from the relativistic limit, it is not clear what the equivalent of ( |2.39| ) should 
be. However, the charge normalization can still be used. Indeed, even for excitations of 
finite energy, since F{—oo\9p) is independent of 9p, the renormalized charge is as before 
^/2g. From ( |3.3|) we conjecture 

/+_(^i,^2) = -iV^expH7rF(^i|^2)][6(^i)]^/'[e(^2)]^/'(-^) ' ' • (3.4) 



Here, e = 4|. This formula agrees with ( p.45|) in the relativistic limit once the change of 



normalization 6{9i — ^2) ^iPi — P2) is taken into account. The last term is fixed by 
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dimensional analysis, together with the charge normalization (it is necessary to off set the 
overall exp[— z7rF(6', 6*)] at coincident rapidities). 

It is interesting to discuss this form factor by getting back to the noise without impu- 
rity. The contribution of the two particle form factor is proportional to (assuming w > 0) 

/ / |/+_(^i,^2)|'5[e(^i)-e(^2)-c^]ciM^2= / —-\f+-ie,e-'[e{e)-u;])\^de. 
Ja J a J a el^/J 

(3.5) 

At first order, we simply set q = ^/2g, 9 = A, to get 



y^gf / = 2guj 

J A 



as requested. At next order, we have to take the other terms into account. Setting 6i 
observe that 

F{e^\e2)=L{A,A){e^-e2) 

=L{A, A){ei-A + A-e2) 
L(A,A) 



eiA) 
-Xuj 



Here we defined 



_ LjA^A) _ LjA^A) 



The contribution of this term to the noise reads 

J A 



\u-e{e) 

2guj I exp 





Xu) In 



X 



Expanding the exponential to first order gives no correction due to In = 0. 

Of course, the form factor contributes a non vanishing term at order uj^. This should be 
offset by the next form-factor involving a pair of voltons and a pair of antivoltons. Indeed, 
while in the relativistic limit with k 7^ 0, such a form factor gave a contribution of order 
uj to the noise, now the residue axiom, reasonably generalized to the non-relativistic case, 
will lead to a form factor similar to ( |2.47| ), but with the term e^*'^'* — 1 replaced by a term 
of order y, measuring the difference between the S-matrix and —1 away from the Fermi 
rapidity. In the computation of the noise, this will give a term of order uj{uj/V)'^. More 
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generally, the process involving n voltons and n antivoltons will give a leading contribution 
of order a;(a;/y)^"'~^. 

Let us now discuss in more details the term uj'^. The integral (|3.5| ) gives for the noise 
a contribution (assuming w > 0) 

S[uj)i^) = 3-^ + C^\ (3.7) 

7T 

where the value of C, will not be needed in what follows. Since the total noise is S{uj) = 
this means that the 4 particle form factor must contribute to the noise by a leading term 

^(u;)W = -Cc^^ (3.8) 

and thus we obtain the sum rule for the four particle form factor 

/ |/++__(^i, ^2, ^3, ^4)|'27r5[e(^i) + e(^2) - e(^3) - e{e^) - u]\\-^ = (3.9) 

The point is, that this sum rule is enough to determine the next term in x"(u;). Indeed, 
the two particle form factor contributes to x" by a term proportional to 

/•oo /•oo 

/ / |/+-(^i, ^2)|'5[e(0i) - e(^2) - ^^]Re[7^(^l)7^*(^2) - l]rf^ici02. (3.10) 

J A J A 

We are now interested in getting the terms of order uj"^ and uj^ in this integral. The term of 
order w"^ was already evaluated in I. For the term of order u;^, since the integral combined 
with the delta function contributes an overall l/uj, there will be two contributions: we can 
either expand the TZ bracket to order and take the leading expression {0{uj'^)) for the 
two particle form-factor, or expand the TZ bracket only to order uj"^ and take the next to 
leading contribution for the two particle form factor {0{uj^)). 

Similarly, the four particle form factor contributes by a term proportional to 

/•oo 

/ |/++__(^i, ^2, ^3, ^4) |'5[e(^i) + e(^2) - e(^3) - <e^) - ^] 

J A 

4 (3.11) 

Re[7^(^l)7^(^2)7^*(^3)7^*(^4) - i]\[dei. 

1=1 

Taking the term of order u)^ in the four particle form factor would necessitate taking the 
order in the TZ bracket, which vanishes. Hence, the only term we need to consider is order 
uj'^ in the TZ bracket, and the leading order (0(0;"*)) in the four particle form factor. Now 
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the point is, that hke in the previous section (paragraph 2.5), the term of order uP' in the 
7^ bracket is actuaUy independent of the rapidities, and simply factors out as a constant. 
The integral that is left is the same as the one appearing in the noise. The same remark 



holds for the two particle form-factor contribution ( 3.10 ), so these two terms cancel out! 
Hence, all what remains to determine the uT' order is the 71. bracket to order cu^, combined 
with the leading expression (0(a;^)) for the two particle form-factor. In other words, the 
expression that was obtained in I is good to get the order: 



-1 Aw uj 

■ ^ ^ Re / dMP'^W{e2)n{e'^)-l\e^^e^'^5{uj-e^^ -e^'^), (3.12) 



Here we have reparametrized e ^ e^, which defines the 9 —* (3 correspondence. Redefining 
n{e) = n{e) we rewrite (|3T2D as 



1 r 

- — 2 — 2-^^ / dx\R{x)7i*{x — uj) — l\dx 

1 

^5— Re / dx[JZ{ujx)TZ* {xuj — uj) — 1]. 

T^^<^ Jo 



(3.13) 



2g7r 



Laborious but straightforward manipulations lead to the expression 



UJ 
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where we have set as usual TZ = e^^ , and primes denote successive derivatives with respect 



+ 
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(3.14) 



to the variable x = e{9). Equation (|3.14|) can be made completely explicit using the exact 
form of the dressed R matrix eq. (5.4) of I, together with the relation between e and 9 
as given in eq. (2.13) of I. No simplification emerges, and at this stage there seems little 
point in giving a more explicit expression of the second term in (|3.14|). 



4. The problem with F = 0, T > 
4.1. Generalities 

When T 7^ 0, it is well known that thermal properties can be computed by evaluating 
correlators in a "thermal ground state", that is, any state that is characterized by the 
equilbrium densities 0], [0. For the problem at hand, restricting once again to g = 1/t, 
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recall that the particles are soliton, antisoliton and the n-breathers n = 1,2, . . .,t — 2. 
Particles will be designated by the generic label j. Introduce the functions solution of 

/oo 
^jk{e-e')\n i + e-^'=(^')/'^ de'. (4.1) 
-00 



k 

Then, equilibrium densities are given by 



and the filling fractions are 

P_3 

The e's diverge at infinity as in a free theory: tj{9) ~ /UjC^, >> T. At minus infinity, 
an interesting consequence of the interaction is that the e's go to a constant of order T. 
Explicitely one has 



e 



(n + 1)2 -l,e'^-''/'^^t- 1,^^-00. (4.4) 



It is well known (see also below) that the e's are excitation energies: destroying (creat- 
ing) a particle of type j decreases (increases) the energy by e^. Hence, due to ([4.4|), there 
are no single particle excitations of arbitrarily low energy: the low energy excitations are, 
as in the voltage case, fully obtained by particle hole pairs (there is a gap of order T for 
any processes that do not conserve the number of particles). Contrarily to the foregoing 
voltage case, such pairs can however be created at any (neighbouring) rapidities since at 
T > 0, the filling fractions are less than one. 

As in the voltage case, we would like to know the form factors of the current. Observe 
that, since only particle hole pairs are involved in the very low energy limit, breathers 
won't appear (since the current operator changes the parity of the number of breathers) 
and we can restrict to soliton/antsoliton excitations. To address the form factors, we need 
as before the renormalized charge and S-matrix. 

4-2. Renormalized S matrix, charge and form-factors 

We start with the shift function. Suppose for instance we add a particle of type k at 
rapidity 9p. The shift of particles of type j obeys, similar to 



/oo 
^jiiO - 0'Me')Fik{6'\6p) = 6,/^k{0) - SjkiO - Op). (4.5) 
-00 
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Observe that, since we consider massless excitations around arbitrary rapidities, we have 
aUowed the additional phase shift St to depend upon 9. 

To proceed, let us assume there is only one type of excitations. Dealing with several 
types makes the notation more cumbersome without changing the argument. Then the 
equation for the shift reads 

/oo 
HO - 9')aO')F{e'\ep) = 5b{d) - 5{0 - Op). (4.6) 
-oo 

The excitation energy for adding this particle is 

/oo 
^e''aO')F{e'\ep). (4.7) 
-oo 

Note that, compared with ( |2.16| ), there is no term arising from —oo since without a voltage. 



the bare (kinetic) energy of particles vanishes in that limit. From ( |4.7| ) we deduce 



dOp J -oo dOp 



Now, from (|4.6|) we have 



'p 

Introduce the kernel, similar to ( p.4| ) 

L{e\e' 

It follows that 



/oo 
- e")^{0")^e" - e') + . 
-oo 

dF{9\9p) _ L{9\9p) _ L{9p\9) 



(4.9) 



d9p aop) m 

and thus 

/oo 
^e^' L{9p\9')d9' 
- oo 



de{9^ 



d9p 



On the other hand, from ( [4.1| ) in the case of a single particle, we deduce 

^ = ^e^ + 1^ HO - (^'n0')%d9\ (4.10) 



so we find 

de de 
d9p d9p 
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hence e{9p) = e{9p) + est, the constant depending on the function St,{9). 

Consider now excitations around some rapidity ^o- Shift the energy and momentum 
such that these excitations have e = p = e{9) — e{9o) = e{9) — e(6'o). In the low energy 
hmit, parametrize e = e^, we are then in a situation similar to excitations around the Fermi 
rapidity in the previous case. Chose the function 5b{9) such that F{9q\9q) = for any 9q. 
Then the analog of k vanishes. To completely determine the form factor of the current, 
we need to determine the charge of excitations. This is simpler than in the case of finite 
voltage. Indeed, since solitons and antisolitons behave in the same way except for their 
charge, we chose the same function Sb for both species, and Fs{—oo,9p) = Fa{—oo,9p). 
So, as many solitons as antisolitons move across the rapidity — oo, resulting in a vanishing 
contribution to the charge. In other words the renormalized charge equals the bare charge. 
Going back to 9 variables, at low energy 9i ~ ^2, 

< 0|j(;^)|^2,^i >-+= -z[e,(^i)]V2[e^(^2)]'/'e^^[^^(^i)-^^(^^)], (4.11) 

where we also used the symmetry (in the absence of applied voltage) between solitons and 
antisolitons eg = 

4.3. DC Conductance in the fractional quantum Hall effect using Kuho formula 
From ( [4.11j ), the u; — limit of the noise without impurity follows 

j dtmim)) 

^^^[e,(^i)] (1 - ^[e,(^2)]) e,(^i)e,(^2)27r5[e,(^i) - -u] + {u^ -u) 



00 
00 



= V..(-oc) ^ l + e-/^l + e-WT + = ^ + ). 

(4.12) 

as required i. 

In the presence of an impurity, we can also compute the noise in a way similar to 
paper I. The LL and LR noise take the same form. The RR noise is a bit different because 
in the (thermal) ground state we have both solitons and antisolitons this time. In fact the 



To compare with the usual form of Nyquist theorem, recall that we have set Tx = 1, so the 

g 
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dimensionless conductance of the pure Luttinger liquid reads G — — 



RR noise is easily seen to be the same as the LL noise here, as is easily proven by changing 
basis \e P\e >| +Q\e Hence 

J dt{{j{t),j{o)}) = ^ Pj{e)[i-mW)\Pm^de. 

Prom this we deduce the formula for the DC conductance using the Kubo formula @ 

1 f°° 

G=7f; UO)[l - UO)¥smm?de, (4.13) 

in agreement with Observe that, should multiple particle excitations be allowed as 
discussed briefly in the voltage case, terms involving higher orders of |Pp would appear in 
G. This seems unlikely in view of the success of formula ( [4.13| ) compared with numerical 
simulations. Also, ( 4.13|) was initially derived using Landauer-Biittiker scattering formula. 



and although it is not entirely clear whether one can apply it to the quasiparticles of 
integrable systems, this alternate derivation still increases our confidence in (|4.13|) . 

4.4- Dynamical susceptibility in the double well problem 

As an application, we would like to consider the low frequency behaviour of the dy- 
namical susceptibility in the dissipative quantum mechanics problem and T ^ 0. We need 
to introduce a dressed reflection matrix in that case, which reads here 



UjiO) = Rj{9)exp 
We define at T ^ 0, 



/oo T 
Jk{e')F,,{9'\9) — lnR,{e') 



(4.14) 



= ^tanh^ / §^e^'''{{S'it),S%0)}). 



(4.15) 



The correlator of the spin anticommutator can be related with the correlator of current 
operators. The argument works exactly like at T = (the propagator with Neumann 
boundary conditions eq. (5.4) of becomes a rational function of trigonometric functions 
at finite T. When a; — 0, it vanishes everywhere except at y = y' , and its integral for y 
running on the interval [0, 1/T] is still equal to 1). Using [0 and the previous form-factors. 



Here, we defined G as 27r times the usual conductance. 
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one finds at low frequency ( using again tlie symmetry between the two species at vanisliing 
bias) 

2T 1 1"^ 

— X"(c^) ^ -^^Re / i[esm (1 - i[eM + ^]) UO) (7^[e(^)]7^*[e(^) dO. 

(4.16) 

Note that ( [4.16|) differs from ( |3.12| ) by factors of 2, due to the fact that excitations involve 
solitons or antisolitons, and then can have positive or negative energy. There is also a factor 
of 2g that differs, due to the fact that the renormalized charge is 1 with a temperature and 
no voltage, while it is ^/2g with a voltage and no temperature. From ( [4.16| ) one finds thus 



0.-0 u.''^ ^ ' 2gV J_„ e.(e) \dO J ^ ' 

In the case g = 1/2, there is no dressing, and one finds 

(^^0 u ^ Jo 1 + 1 + e ^/^ \x'^ + T^J 

so, after using the transformation 
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X + 1 X — 1 x^ — 1 
one finds, using standard integral representations of the t/j function. 



hm — X {^) 



in agreement with the result of [ p!5| . 



Tb V27rT 2 J 2^ \27rT 2 



(4.19) 



As small temperature, the sums are dominated by particles whose rapidities approach 
— oo. In that limit, one checks that ^ InT?. ~ f^P? from which we recover x"{^) ~ gi-KTsY 
at T = 0, as shown in I. 

Eq. ( [l.lTp can be made as explicit as necessary. Using the relation 



one finds 



in7e,(^) = -ini?,(^) + J] j L,,ie\e')—inR,ie'). 



Here, one has 

/oo 
-oo 

Since functions ^ are easy to obtain by numerical solution of the thermodynamic Bethe 
ansatz equations, it is only a technical matter to determine ( [1.17|) . 
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5. Conclusion 



We hope that the computations in this paper and the previous one open a possible 
direction for the study of correlators at finite T and V in quantum impurity problems. Our 
results ( |3.14| ),( ^T7| ), although completely explicit, are unfortunately quite cumbersome: 
pushing the method further and getting results at arbitrarily large frequency looks much 
more involved than in the T = case, V = case (the situation looks more favorable for 
problems which are massive in the bulk, especially at low temperature 0). Nevertheless, 
interesting questions are already raised by this approach: in particular, the results of the 
present paper and those found in I do not agree with those of perturbative methods. It 
would be very useful to investigate the matter further in experiments or numerical studies 
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